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HIGHLIGHTS 


•  Stress— strain  relationship  was  measured  to  determine  thickness  and  porosity  of  DM. 

•  Thermal  conductivity  was  found  highly  dependent  on  compression  and  water  content. 

•  Theoretical  prediction  of  the  maximum  thermal  conductivity  was  developed. 

•  Internal  cell  temperature  difference  was  estimated  as  a  function  of  saturation. 

•  PCI  flow  is  more  influential  on  multiphase  transport  with  DM  at  lower  saturation. 
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An  experimental  study  to  investigate  the  through-plane  thermal  conductivity  of  three  different  diffusion 
media  (DM)  used  in  polymer  electrolyte  fuel  cells  (PEFCs)  as  a  function  of  compression  (from  0.1  MPa  to 
2  MPa)  and  saturation  (from  0  to  25%)  was  performed.  Additionally,  measurements  to  determine  the 
stress— strain  relationship  for  the  materials  were  made  using  an  optical  microscope.  Both  compression 
and  water  content  had  a  significant  impact  on  the  through-plane  thermal  conductivity,  which  should  be 
accounted  for  in  multiphase  modeling  efforts.  An  analytical  expression  for  the  theoretical  maximum  of 
the  through-plane  thermal  conductivity,  as  a  function  of  both  compression  and  saturation,  was  devel¬ 
oped  to  help  understand  the  nature  of  liquid  connectivity  in  saturated  pores.  Additionally,  a  relationship 
was  developed  to  predict  actual  thermal  conductivity  of  the  tested  materials  as  a  function  of  both 
compression  and  saturation  based  on  experimentally  measured  data. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Carbon  paper  is  used  as  the  diffusion  medium  (DM)  in  polymer 
electrolyte  fuel  cells  (PEFCs)  due  to  its  suitable  permeability,  elec¬ 
tronic  conductivity,  thermal  conductivity,  and  mechanical  strength. 
Thermal  conductivity  is  a  particularly  important  parameter,  due  to 
the  interplay  between  heat  and  water  management  [1—4]  and  the 
perceived  dominance  of  conduction  heat  transfer  in  PEFCs.  Water 
generated  in  PEFCs  that  condenses  into  liquid  phase  must  be 
removed  to  prevent  loss  of  performance  during  operation  and  po¬ 
tential  freeze-thaw  degradation  [5-8].  Polytetrafluoroethylene 
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(PTFE)  is  normally  used  to  improve  hydrophobicity  of  the  DM  and 
reduce  liquid  accumulation.  A  highly  hydrophobic  microporous 
layer  (MPL)  is  often  utilized  between  the  DM  and  catalyst  layer  (CL). 
The  MPL  is  of  great  importance  to  mitigate  flooding,  to  decrease 
electrical  contact  resistance,  protect  the  membrane  from  damage 
from  DM  fibers,  and  prevent  dryout  from  excessive  vapor  removal. 

Various  experimental  approaches  have  been  employed  to 
measure  the  thermal  conductivity  of  unsaturated  DM  [9-14]. 
Khandelwal  and  Mench  [11  measured  carbon  paper  and  Nation® 
membranes  for  PEFCs.  They  reported  that  Sigracet®  20  wt.%  PTFE 
treated  carbon  paper  had  a  thermal  conductivity  of 
0.22  ±  0.04  W  m-1  K_1  and  Toray  paper  had  a  thermal  conductivity 
of  1.80  it  0.27  W  nrT1  K_1.  Burheim  et  al.  [12]  reported  an  increase 
of  thermal  conductivity  of  artificially  aged  SGL  gas  diffusion  media 
while  PTFE  content  decreases.  They  also  measured  SolviCore  DM  at 
different  compression  pressures.  They  reported  that  thermal 


G.  Xu  et  al.  /  Journal  of  Power  Sources  256  (2014)  212-219 


213 


conductivity  increases  as  compression  pressures  grow  and  an  in¬ 
crease  in  thermal  conductivity  when  the  DM  is  saturated  [13]. 
Zamel  et  al.  [9]  measured  the  through-plane  thermal  conductivity 
of  unsaturated  Toray  paper  in  a  temperature  range  from  -50  °C  to 
120  °C.  A  thermal  conductivity  of  0.8-1.8  W  m-1 1<-1  was  reported 
at  high  compression,  and  0.2-0.4  W  m"1  K  1  at  low  compression. 
Although  measurement  of  dry  media  is  useful,  in  operation,  PEFCs 
commonly  have  DM  saturation  at  levels  up  to  30%  [15-17]. 
Therefore,  knowledge  of  partially  saturated  conductivity  is 
important. 

Several  researchers  also  proposed  numerical  models  to  calculate 
thermal  conductivity  of  DM  [18-20].  In  order  to  better  understand 
the  mechanism  of  water  transport  and  its  effect  on  performance 
and  durability,  Bazylak  provided  an  overview  of  recent  de¬ 
velopments  in  liquid  water  visualization  with  PEFCs  [21  .  Yablecki 
et  al.  [22]  reported  an  increase  of  20.8%  in  through-plane  thermal 
conductivity  for  their  gas  diffusion  layer  (GDL)  modeling  domain 
with  24.4%  saturation.  However,  measurement  of  the  thermal 
conductivity  of  diffusion  media  as  a  function  of  both  compression 
and  a  full  range  of  saturation  has  not  yet  been  fully  developed,  so 
that  the  impact  on  transport  can  be  fully  understood,  which  is  the 
motivation  for  this  study.  In  this  work,  measurements  of  stress- 
strain  relationship  and  thermal  conductivity  as  a  function  of  satu¬ 
ration  were  performed.  A  theoretical  analysis  of  the  maximum 
thermal  conductivity  as  a  function  of  saturation  has  been  devel¬ 
oped  to  glean  insight  into  the  internal  distribution  and  connectivity 
of  liquid  in  the  media.  Water  transport  within  the  saturated  DM 
was  also  investigated  via  calculations  of  phase-changed-induced 
(PCI)  flow  and  its  impact  as  a  function  of  saturation  [23].  The  re¬ 
sults  of  this  work  should  be  useful  to  understand  and  more  pre¬ 
cisely  model  thermal  transport  in  operating  PEFCs  with  realistic 
saturation  distributions. 


2.  Experimental  setup 

Diffusion  media  were  uniformly  compressed  between  two 
25  mm  thick  stainless  steel  cylinders  with  known  thermal  con¬ 
ductivity.  This  design  enabled  determination  of  the  stress-strain 
relationship  and  the  through-plane  thermal  conductivity.  Sche¬ 
matic  diagrams  of  the  experimental  apparatus  are  shown  in  Figs.  1 
and  2.  A  load  cell  was  used  to  measure  compression  pressure  on  the 
diffusion  media.  All  tests  were  performed  at  room  temperature, 
which  varied  from  21  °C  to  27  °C  throughout  the  course  of  testing. 


3.  Design  and  instrumentation 

A  Hot  Disk  TPS2500S  Thermal  Constants  Analyser  (Hot  Disk  AB, 
Gothenburg,  Sweden  &  ThermTest  Inc.,  Fredericton,  Canada)  was 
used  to  measure  the  thermal  conductivity  of  all  materials  tested  in 
this  work.  Based  on  the  theory  of  the  Hot  Disk  Transient  Plane 
Source  (TPS)  technique,  a  TPS  sensor  in  the  shape  of  a  double 
spiral  was  utilized  to  take  measurements  within  ten  seconds, 
effectively  decreasing  unwanted  water  vaporization  during  the 
tests.  The  TPS  sensor  acted  both  as  a  heat  source  for  increasing  the 
temperature  of  the  samples  and  a  resistance  thermometer  for 
recording  the  time  dependent  temperature  increases  [24].  The 
reproducibility  of  measured  thermal  conductivity  using  this 
technique  was  observed  to  be  ±3%.  A  TPS  7280  thin  film  sensor 
with  a  sensor  diameter  of  29.34  mm  was  used  for  all  testing.  This 
large  format  sensor  provides  a  uniform  heat  flux  across  the 
measurement  area. 

4.  Measurement  and  uncertainty  analysis 

Round,  53  mm  diameter  samples  were  used  for  stress-strain 
measurement.  Mitsubishi  Rayon  Corp.  Grafil  U-105  (MRC  105), 
SGL  Sigracet®  25  BC,  and  General  Motors  (GM)  Experimental  virgin 
samples  were  used  in  this  study.  GM  Experimental  samples  are 
materials  that  satisfy  the  2015  US  DOE  performance  targets  while 
enabling  progress  toward  the  automotive  fuel  cell  system  cost 
target  of  $30/kWe  as  developed  under  DOE  Award  Number  DE- 
EE0000470.  To  reduce  measurement  error,  two  samples  of  MRC  105, 
SGL  25  BC,  or  GM  Experimental  were  layered  for  each  stress-strain 
measurement.  Images  of  compressed  samples  were  captured  using 
a  zoom  stereomicroscope  (Unitron  ZST)  and  digital  CCD  camera 
(Olympus  DP  73)  with  optical  resolution  of  1600  x  1200  pixels. 
Deflections  of  the  DM  in  the  images  were  measured  using  digital 
imaging  software  cellSens  Standard. 

In  this  study,  the  DM  stress-strain  relationship  was  optically 
measured,  as  shown  in  Fig.  3.  The  thickness  versus  compression 
relationship  was  then  used  to  calculate  the  thickness  used  in  the 
determination  of  compressed  DM  thermal  conductivity.  Porosity  is 
of  great  importance  to  calculate  the  actual  saturation  of  the  porous 
media.  As  compression  increases,  the  void  volume  in  DM  is 
reduced,  decreasing  effective  porosity  according  to  the  following 
equation: 


Digital  CCD  Camera 
Olympus  DP  73 


1600x1200 
lxl,  3CCD  mode 


Fig.  1.  Schematic  diagram  of  the  experimental  apparatus  used  to  optically  measure 
stress-strain  in  thin  film  materials  (not  to  scale). 


Fig.  2.  Schematic  diagram  of  the  experimental  apparatus  used  to  measure  thermal 
conductivity  of  thin  film  materials  (not  to  scale). 
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Fig.  3.  Micrograph  of  a  DM  under  compression. 


where  0*  is  the  compressed  porosity  of  DM,  0  is  the  initial  porosity 
of  DM,  and  e  is  the  total  strain  on  the  compressed  DM.  It  should  be 
noted  that  during  compression  to  2  MPa,  the  MPL  is  essentially 
incompressible  when  layered  as  a  composite  with  the  much  less 
stiff  diffusion  media,  so  that  the  load  deformation  can  be  consid¬ 
ered  on  the  DM  only.  Table  1  shows  initial  porosity  and  compressed 
porosity  of  three  different  diffusion  media  at  compression  of  2  MPa. 


4  A.  Saturation  of  the  DM 

Unsaturated  weight  and  thickness  were  measured  prior  to  the 
saturation  process  using  a  Sartorius  TE214S  analytical  balance.  In 
this  study,  porous  diffusion  media  were  artificially  saturated  by 
immersion  in  de-ionized  (DI)  water  in  a  sealed  container  while  a 
vacuum  was  applied.  A  sonic  bath  was  utilized  to  remove  residual 
bubbles  during  immersion  and  saturation  [25]. 

The  thickness  at  compression  of  2  MPa  was  calculated  based  on 
the  stress  strain  relationship  and  the  following  equation: 


where  l  is  the  thickness  of  compressed  DM,  is  the  thickness  of 
measured  DM,  e  is  the  total  strain  on  compressed  DM,  e°  is  the 
strain  on  measured  DM. 

The  saturated  weight  of  the  diffusion  medium  was  calculated 
prior  to  saturation  based  on  the  following  equation: 

fftwet  —  tn^j-y  +  spwa|-erA/j  0  (3) 

where  mwe t  and  m^y  are  weights  of  saturated  and  unsaturated 
diffusion  media,  respectively,  s  is  the  saturation  of  DM,  pWater  is  the 
density  of  water,  A  is  the  area  of  DM,  l\  is  the  thickness  of  the 
compressed  macroporous  layer,  0*  is  the  compressed  porosity  of 
the  DM. 

For  MRC  105  DM,  without  MPL,  the  thickness  of  the  DM  l\  is  the 
thickness  of  compressed  DM.  For  GM  Experimental  DM  shown  in 
Fig.  4,  with  thin  MPL,  the  thickness  of  the  compressed  macroporous 


Table  1 

Calculated  change  in  porosity  at  compression  of  2  MPa. 


DM 

Initial  porosity 

Compressed  porosity 

GM  Exp. 

82.0% 

67.8%  ±  1.2% 

SGL  25  BC 

80.0% 

64.8%  ±  0.7% 

MRC  105 

83.2% 

79.1%  ±  0.2% 

layer  was  almost  the  same  as  the  thickness  of  compressed  DM  since 
the  MPL  is  relatively  thin  compared  to  the  macroporous  layer.  For 
SGL  25  BC,  with  a  thicker  MPL,  the  thickness  of  the  compressed 
macroporous  layer  was  used  to  calculate  saturation  instead  of  the 
thickness  of  compressed  DM  since  the  highly  hydrophobic  MPL  was 
assumed  to  be  unsaturated.  The  compressed  porosity  of  DM  was 
calculated  as  shown  in  Table  1. 

Since  some  water  can  be  squeezed  out  of  the  diffusion  media 
during  compression  loading,  the  average  decrease  in  saturation 
level  during  compression  loading  was  calculated  and  compensated 
prior  to  saturation  process.  After  each  test  series,  an  updated 
saturated  weight  was  immediately  measured.  Saturation  was 
calculated  based  on  the  following  equation: 

s  =  mwet  -  mdy  (4) 

P  water'll 0 


4.2.  Unsaturated  and  saturated  DM  thermal  conductivity 
measurement 

Thermal  conductivity  was  measured  with  the  through-plane 
test  module.  Two  specimens  were  required  for  each  test.  As 
shown  in  Fig.  2,  one  sample  was  placed  above  the  Hot  Disk  sensor. 
The  other  was  placed  underneath  the  sensor.  The  assembly  was 
firmly  clamped  between  two  25  mm  thick  stainless  steel  cylinders 
with  known  high  thermal  conductivity.  Samples  were  consecu¬ 
tively  tested  at  compressions  of  0.1  MPa,  0.5  MPa,  1  MPa,  1.5  MPa 
and  2  MPa.  Thirty  minutes  of  pre-test  time  was  taken  before  each 
test  to  ensure  steady  state  conditions.  The  actual  measurement 
time  was  10  s  at  a  heating  power  of  2.2  W.  The  ambient  temper¬ 
ature  in  this  study  varied  from  21  °C  to  27  °C,  which  should  not 
make  a  significant  impact  on  results  based  on  previous  studies 

ini. 

In  the  experiments  for  determination  of  saturated  DM  thermal 
conductivity,  the  compression  pressure  was  set  at  2  MPa,  which  is 
within  the  normal  range  for  a  fuel  cell  assembly.  Samples  were 
initially  saturated  to  a  high  saturation  level  and  allowed  to  evap- 
oratively  dry  until  the  target  saturation  level  was  reached. 


4.3.  Measurement  validation 

The  thickness  measurement  was  validated  with  a  stage  micro¬ 
meter  for  transmitted  light  (Ted  Pella,  Inc.).  The  measured  thickness 
was  1.89%  lower  than  the  actual  value,  which  is  within  the  expected 
experimental  uncertainty. 

Thermal  measurement  was  validated  with  a  thin  disk  of  Teflon® 
PTFE,  SGL  25  BC  DM  and  stainless  steel.  The  measured  PTFE  thermal 
conductivity  was  determined  to  be  0.240  ±  0.025  W  m  1  K-1 
compared  to  the  Teflon®  PTFE  standard  thermal  conductivity  of 
0.259  ±  0.006  W  m-1  K-1.  The  measured  SGL  25  BC  thermal  con¬ 
ductivity  was  0.214  ±  0.004  W  m-1  K_1  at  compression  of  0.5  MPa 
and  0.356  ±  0.018  W  m-1  I<-1  at  compression  of  1  MPa.  Ramousse 
et  al.  [26]  reported  a  thermal  conductivity  of  0.26-0.34  W  m-1 K_1. 
Khandelwal  and  Mench  [11]  reported  thermal  conductivity  of 
0.31  ±  0.06  W  m-1  K-1.  The  measured  stainless  steel  thermal 
conductivity  was  13.380  ±  0.036  W  m-1  K-1.  Hot  Disk  AB,  a  Sweden 
company  which  develops,  designs,  manufactures,  sells,  installs  and 
services  Thermal  Constants  Analyzer,  reported  an  average  thermal 
conductivity  of  13.7418  W  m-1  K^1  with  less  than  0.15%  standard 
deviation  based  on  these  results,  the  test  system  was  deemed  to  be 
within  acceptable  limits  of  accuracy  [27]. 
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Fig.  4.  Configuration  of  DM  a)  GM  Exp.  MPL  coated  on  GDL  b)  SGL  25  BC  MPL  coated  on  GDL  c)  MRC  105  GDL  without  MPL  coating  d)  GM  Exp.  face  of  macroporous  layer  e)  SGL  25 
BC  face  of  macroporous  layer  f)  MRC  105  face  of  macroporous  layer  g)  GM  Exp.  cross  section  h)  SGL  25  BC  cross  section  i)  MRC  105  cross  section. 


5.  Results  and  discussion 

5.1.  SEM  images 

The  characteristic  structures  of  the  cross-section,  the  face  of  the 
macroporous  layers  and  the  face  of  the  microporous  layers  are 
shown  in  Fig.  4  using  scanning  electron  micrographs  (HITACHI 
TM3000). 

The  micrographs  show  that  SGL  25  BC,  MRC  105  and  GM 
Experimental  DM  primarily  have  planar  structures  of  fibers  aligned 
in  multiple  layers.  The  presence  of  PTFE  can  be  observed  between 
the  fibers.  SGL  25  BC  has  a  relatively  thick  MPL  (46%  of  DM 


Stress  (MPa) 

Fig.  5.  Experimental  data  and  curve  fit  lines  for  measured  stress-strain  relationships. 


thickness),  which  is  essentially  incompressible  under  the  loads  in 
testing.  The  GM  Experimental  samples  have  a  much  thinner  MPL, 
and  the  MRC  105  samples  do  not  have  a  MPL.  In  the  GM  Experi¬ 
mental  DM,  a  carbon-fiber  skeleton  filled  with  carbon  powder  is 
observed,  as  shown  in  Fig.  4. 

5.2.  Stress-strain  relationship 

Fig.  5  presents  three  fitting  curves  of  the  measured  stress-strain 
deflection  for  GM  Experimental,  SGL  25  BC  and  MRC  105.  The 
stress— strain  relationship  for  the  DM  is  essential  to  adjust  for 
thermal  conductivity  measurements  for  actual  thickness.  For  GM 
Experimental  DM,  the  strain  sharply  increases  15%  during  initial 
compression  loading  (0.1  MPa)  and  nonlinearly  approaches  44%  at 
compression  of  2  MPa.  The  effective  compressed  porosity  is 
reduced  from  82.0%  to  67.8%  during  compression  to  2  MPa.  For  SGL 
25  BC  MPL,  the  total  strain  is  calculated  to  be  4.9%.  For  SGL  25  BC 
macroporous  layer,  the  strain  reaches  10%  during  the  initial 
compression  loading,  and  approaches  31%  resulting  in  reduced 
porosity  from  80.0%  to  64.8%.  MRC  105  samples  do  not  have  as 
much  strain  during  the  initial  compression  loading.  Following  a  9% 
strain  at  compression  of  0.25  MPa,  MRC  105  samples  deflect  linearly 
up  to  19.5%  at  compression  of  2  MPa.  A  nonlinear  curve  fitting  was 
performed  using  an  exponential  association  model: 


Table  2 

Measured  stress-strain  curve  fit  functions. 


DM 

Experimental  stress-strain  curve  fit  functions 

GM  Exp. 

e  =  0.14429(1  -  e"^0  04031) 

+  0.31607(1  -  e_<7/0-70534) 

SGL  25  BC 

e  =  0.29558(1  -  e"^1-68126) 

+  0.10073(1  -  e-ff/0  02651) 

MRC  105 

e  =  0.07055(1  -  e~al0A2226) 

+  1.79812  *  108(1  -  e(7/2-98335*109) 
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Fig.  6.  Measured  thermal  conductivity  of  unsaturated  DM  as  a  function  of 
compression. 


Y ;•=  a-i  (l  -  e-£)  +  a2  (l  -  e“£)  (5) 

Table  2  and  Fig.  5  show  the  three  experimental  stress-strain 
curve  fitting  models  of  GM  Experimental,  SGL  25  BC  and  MRC  105, 
respectively. 

A  residual  strain  can  be  observed  when  the  compressive  load  is 
released  due  to  plastic  deformation.  In  this  study,  the  focuses  are 
thickness  measurement  and  prediction,  and  thermal  conductivity 
measurements  were  only  performed  on  virgin  materials  under  first 
compression.  Therefore,  successive  compress-release  cycles  were 
not  performed. 

5.3.  Thermal  conductivity  measurement  of  unsaturated  samples 

Fig.  6  and  Table  3  show  measured  thermal  conductivity  aver¬ 
aged  across  five  tests  using  material  from  the  same  lot  to  reduce 
manufacturing  variability  and  linear  fitting  functions  derived  from 
unsaturated  DM  under  various  compressions  up  to  2  MPa.  Contact 
resistance  between  the  samples  and  the  sensor  has  been  removed 
from  the  measured  data  by  the  device  as  the  influence  of  contact 
resistance  only  occurs  during  the  initial  portion  of  the  transient 
measurement  [28,29].  Thermal  conductivity  was  relatively  lower 
when  low  compression  pressure  was  applied  because  of  the  pres¬ 
ence  of  the  internal  contact  resistance  and  high  porosity  in  the  DM. 
With  increased  compression,  internal  contact  resistance  and 
porosity  decreases,  and  thermal  conductivity  increases.  The  ther¬ 
mal  conductivity  of  unsaturated  GM  Experimental  DM  mono- 
tonically  increases  with  respect  to  compression  pressure  from 
0.1  MPa  to  2  MPa  by  440%.  The  measured  thermal  conductivity 
increased  about  80%  when  compression  pressure  varied  from 
1.0  MPa  to  1.5  MPa.  Since  the  strain  almost  reached  a  maximum, 
carbon  powder  and  carbon  fiber  were  more  connected  with  each 
other  reducing  internal  conductive  resistance  and  forming 
enhanced  heat  transfer  pathways.  The  thermal  conductivity  of 


Table  3 

Measured  unsaturated  thermal  conductivity  fit  functions. 


DM 

Experimental  unsaturated  thermal 
conductivity  linear  fit  functions 

GM  Exp. 

/<dry  =  0.136  +  0.390(7 

SGL  25  BC 

/<dry  =  0.07121  +  0.31777(7 

MRC  105 

kdTy  =  0.13209  +  0.02349(7 

unsaturated  SGL  25  BC  and  MRC  105  DM  linearly  increases  with 
respect  to  compression  loading. 

5.4.  Thermal  conductivity  measurement  of  saturated  samples 

Fig.  7  and  Table  4  presents  measured  thermal  conductivity  and 
fit  functions  of  saturated  GM  Experimental,  SGL  25  BC  and  MRC  105 
DM  at  compression  of  2  MPa.  As  expected,  water  content  has  a 
significant  effect  on  the  increment  of  thermal  conductivity  of  DM 
since  it  replaces  air.  The  thermal  conductivity  of  MRC  105  sharply 
increases  nearly  62%  from  0  to  6%  saturation  and  gradually  ap¬ 
proaches  theoretical  maximum  as  saturation  levels  increase.  The 
theoretical  maximum  shown  is  developed  from  a  simple  analytical 
model  that  was  derived  for  comparative  purposes  only,  and  is 
described  later  in  this  work.  As  shown  in  Fig.  7,  water  content  ob¬ 
tained  in  GM  Experimental  samples  has  an  obvious  influence  on 
thermal  conduction  with  saturation  from  0  to  10%.  Similar  to  MRC 
105,  the  thermal  conductivity  of  SGL  25  BC  gradually  increases  35% 
as  saturation  level  increases  from  0  to  25%.  As  saturation  increases, 
the  measured  thermal  conductivity  approaches  the  theoretical 
maximum  thermal  conductivity  as  a  function  of  saturation  since 
more  connectivity  of  solid  and  liquid  components  occurs  in  the 
compressed  DM. 

5.5.  Prediction  of  thermal  conductivity  of  saturated  samples  based 
on  experimentally  measured  results 

Figs.  6  and  7,  Tables  3  and  4  show  a  good  agreement  on  the 
thermal  conductivity  of  unsaturated  DM  at  compression  of  2  MPa. 
The  intercepts  of  saturated  thermal  conductivity  curve  fit  functions 
in  Table  4  represents  thermal  conductivity  of  unsaturated  DM  at 
compression  of  2  MPa,  which  can  be  substituted  by  the  unsaturated 
thermal  conductivity  linear  fit  functions  in  Table  3,  yielding  the 
prediction  of  the  thermal  conductivity  of  saturated  DM  as  a  func¬ 
tion  of  both  compression  and  saturation,  which  is  presented  in 
Table  5. 

5.6.  Theoretical  maximum  prediction  of  thermal  conductivity  using 
parallel  model  and  series-parallel  model 

It  is  well  known  that  DM  used  in  PEFC  is  a  heterogeneous  ma¬ 
terial  with  carbon  fibers  or  cloth  and  heat  conduction  is  a  dominant 


Max  &  Min  GM  Exp 


Fig.  7.  Measured  and  predicted  maximum  and  minimum  thermal  conductivity  as  a 
function  of  saturation  for  various  DM  at  2  MPa  compression. 
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Table  4 

Measured  saturated  thermal  conductivity  fit  functions. 


DM 

Saturated  thermal  conductivity  curve  fit  functions 

GM  Exp. 
SGL  25  BC 

MRC  105 

J<wet  =  0.9071  +  0.44496(1  -  e~sl008874) 
kwet  =  0.70051  +  1.44259s 

fcwet  =  0.1872  +  0.52403(1  -  e~s/a6967)  +  0.11151(1  -  e-s/003477) 

Table  5 

Predicted  thermal  conductivity  of  saturated  DM  as  a  function  of  compression  and 
saturation  based  on  experimental  measurement. 


DM 

Predicted  thermal  conductivity  of  saturated  DM 

GM  Exp. 

SGL  25  BC 

MRC  105 

kwet  =  0.136  +  0.390(7  +  0.44496(1  -  e-s/008874) 

kwet  =  0.07121  +  0.3 1 777(7+1 .44259s 

kwet  =  0.13209  +  0.02349(7  +  0.52403(1  -  e^s/a6967) 

+  0.11151(1   e-s/a°3477) 

heat  transfer  mode  through  the  DM.  Many  models  have  been 
described  in  literature  18,19,30-33]  to  predict  the  thermal  con¬ 
ductivity  of  porous  diffusion  media.  But  saturation  was  not 
considered  in  those  models,  primarily  due  to  lack  of  experimental 
data  now  available. 

A  simplified  model  meant  only  to  examine  limiting  behavior 
was  developed  in  this  study  with  basic  components  (carbon,  water 
and  air)  connected,  and  parallel  to  the  heat  flow,  as  shown  in  Fig.  8. 
This  represents  the  theoretical  maximum  of  thermal  conductivity 
of  porous  diffusion  media  at  compression  of  2  MPa,  if  all  saturation 
was  connected  and  aligned. 


k\  —  0  (kw  —  ka)s  ^1  —  0  ^  kc  +  0  ka 

(6) 

II 

(7) 

k-  1 

h  ,  bli 

ki  +  k2 

(8) 

For  SGL  25  BC: 

4.12232s +  8.22939 

(9) 

s  + 8.68781 

For  MRC  105: 


k  =  0.4398s  +  0.437  (10) 

For  GM  Experimental: 

k  =  0.376968s +  1.29783  (11) 


Fig.  8.  Illustration  of  simple  model  to  predict  theoretical  maximum  thermal  conduc¬ 
tivity  of  partially  saturated  porous  diffusion  media. 


where  kw  is  the  thermal  conductivity  of  water,  kc  is  the  thermal 
conductivity  of  carbon,  ka  is  the  thermal  conductivity  of  air,  l\  is  the 
thickness  of  compressed  macroporous  layer,  and  l  is  the  thickness 
of  compressed  DM. 

The  theoretical  maximum  and  minimum  (i.e.  measured  dry) 
prediction  of  thermal  conductivity  using  this  limiting  model  is 
presented  in  Fig.  7.  A  significant  difference  between  measured 
unsaturated  thermal  conductivity  and  theoretical  maximum  was 
observed,  since  the  planar  structures  of  the  fibers  in  all  of  three  DM 
are  orthogonal  to  the  direction  of  heat  flow.  At  low  saturation, 
water  droplets  begin  to  connect  along  multiple  fibers  among 
different  sub-layers  to  create  connected  heat  conduction  pathways. 
These  pathways  significantly  increase  thermal  conductivity.  As 
saturation  increases,  multiple  pores  fill  with  water  forming  more 
connected  heat  conduction  pathways,  approaching  theoretical 
limits.  The  theoretical  limit  represents  a  condition  in  which  all 
water  in  the  diffusion  medium  is  aligned  and  connected  in  the 
direction  of  heat  flow,  maximizing  conductivity  for  a  given  satu¬ 
ration.  The  fact  that  the  measured  curve  would  fail  to  completely 
merge  with  the  maximum  is  expected,  since  a  certain  degree  of 
tortuosity  in  the  DM  exists,  and  the  water  will  not  be  completely 
connected. 


5.7.  Theoretical  maximum  thermal  conductivity  as  a  function  of 
both  compression  and  saturation 

The  compressed  porosity  is  a  function  of  strain,  and  strain  is  a 
function  of  compression,  therefore  compressed  porosity  can  be 
described  as  a  function  of  compression. 

Substituting  Equation  (1)  in  Equation  (6)  yields: 

k\  =  ^1  —  - — [(/<w  —  ka)s  —  kc  +  /<a]  +  kc  (12) 

For  SGL  25  BC  material: 

/  i  * 

(f>  =  0.800,  4  =  -++1  -46%)  (13) 

e*  =  0.1427(l  -e~&)  +  1.88643*1013(l  -  e~ (14) 
For  MRC  105  material: 

0  =  0.832,  h  =1  e  =  e  (15) 

For  GM  Experimental  material: 

0  =  0.820,  l\  ~  /,  e  **  £  (16) 

where  £*  is  the  strain  on  the  macroporous  layer,  0  is  the  initial 
porosity  of  DM. 

Substituting  Equations  (7),  (8),  (13)  and  (14)  and  Table  6  in 
Equation  (12),  the  theoretical  prediction  of  maximum  thermal 
conductivity  for  SGL  25  BC  was  solved: 


Table  6 

Material  properties  of  DM  used  in  this  study. 


DM 

<£* 

(%) 

kw 

(Witt1  K1) 

kc 

(Wm  1  K-1) 

/<a 

(Wm  1  K1) 

hi' 

GM  Exp. 

67.8 

0.58 

3.98 

0.024 

=  1 

MRC  105 

79.1 

0.58 

1.7 

0.024 

1 

SGL  25  BC 

64.8 

0.58 

1.7 

0.024 

0.412 

^Compressed  to  2  MPa. 
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Fig.  9.  Estimated  temperature  difference  across  a  cathode  DM  in  a  fuel  cell  operating 
at  0.6  V,  1.0  A  cm"2  and  353  K. 


heat  flux  calculations,  however,  variation  in  saturation  (thus  cell 
temperature  gradients)  will  also  influence  water  management  via 
PCI  flow. 

Since  PCI  flow  is  governed  by  cell  temperature  gradients,  this 
relationship  was  then  correlated  to  a  saturation  pressure  difference 
[36].  A  PCI  molar  flow  rate  can  be  calculated  as  a  function  of 
saturation  via  an  effective  vapor  diffusion  coefficient  DJff  and  this 
saturation  pressure  difference  [23], 

Df  =  D„0  m  -  s)]1-5  (-Q  "  (21 ) 

where  DVo  is  the  diffusivity  at  absolute  temperature  To  and  pressure 
p0  with  a  constant  n  of  1.81  [37].  The  diffusivity  DVo  was  taken  to  be 
0.226  *  1(T4  m2  s'1  at  273.15  Kand  101  kPa.The  PCI  molar  flux/ was 
then  determined  using  Fields  law  for  diffusion 

J  =  Df g  (22) 


k  _  85(e  —  1)  (695e*S  +  30e*  —  556s  —  449) 

-  w  ^ 

w  =  106250e2  +  34750ee*S  -  162125ee*  -  27800eS  -  14000e 

-  1 8765  (e*) 2s  +  1 1 3940 (/) 2  +  1 7792 £*s  -  1 1 0980e* 

-  2224s  +  55575 

(18) 

Substituting  Equations  (7),  (8)  and  (15)  and  Table  6  in  Equation 
(12),  the  theoretical  prediction  of  maximum  thermal  conductivity 
for  MRC  105  was  solved: 

e  =  0.07055(l  —  e-™®)  +  1.79812*10s(l  - 

(19) 

Substituting  Equations  (7),  (8)  and  (16)  and  Table  6  in  Equation 
(12),  the  theoretical  prediction  of  maximum  of  thermal  conduc¬ 
tivity  for  GM  Experimental  was  solved: 

.  0.556cs  +  0.024c  -  0.462592s  -  0.73608 


5.8.  The  effect  of  thermal  conductivity  on  multiphase  mass 
transport 

It  has  been  shown  experimentally  as  well  as  numerically  that  a 
temperature  gradient  across  diffusion  media  induces  additional 
water  transport  due  to  phase  change,  this  phase-change-induced 
(PCI)  flow  has  also  been  referred  to  as  the  heat  pipe  effect 
[23,34,35].  These  effects  are  significant  in  terms  of  the  overall  water 
balance  of  a  cell.  In  order  to  assess  the  influence  of  PCI  flow,  the 
developed  semi-empirical  thermal  conductivity  relationships  were 
used  to  calculate  a  ratio  between  PCI  flow  and  water  generated 
from  the  electrochemical  reaction.  First,  a  temperature  difference 
across  each  DM  sample  was  calculated  as  a  function  of  saturation. 
Using  the  cathode  catalyst  layer  heat  generation  term  and  Fourier’s 
law  of  conduction,  for  a  single  cell  operating  at  0.6  V,  1.0  A  cm-2, 
353  K,  101  kPa  back  pressure,  the  temperature  difference  across 
each  DM  was  estimated  as  a  function  of  saturation.  This  result, 
shown  in  Fig.  9,  indicates  that  saturation  must  be  considered  to 
ensure  computational  accuracy.  Evidently  this  holds  true  for  cell 


c  =  ^  (23) 

where  C  is  the  local  water  concentration,  which  is  a  function  of 
the  saturation  pressure  Psat,  ideal  gas  constant  R,  and  tempera¬ 
ture  T.  Under  the  chosen  cell  operating  conditions,  the  calculated 
ratio  of  PCI  molar  flow  to  the  molar  flow  of  water  generated  by 
reaction  was  then  determined  as  a  function  of  saturation,  as  seen 
in  Fig.  10. 

The  significance  of  this  result  is  that,  at  lower  saturation,  PCI 
flow  is  more  influential  on  multiphase  transport  within  fuel  cell 
DM.  The  results  also  indicate  that  thermal  resistance  plays  an 
integral  role  in  the  overall  water  balance  and  transport  in  the 
cell,  which  has  already  been  demonstrated  in  the  literature. 
When  this  ratio  is  near  unity,  it  indicates  that  PCI  flow  is  capable 
of  transferring  all  water  generated  within  the  catalyst  layer 
across  the  cathode  DM,  which  agrees  with  the  experimental  re¬ 
sults  of  Owejan  et  al.  [38].  Thus,  in  order  to  properly  account  for 
these  transport  effects  on  overall  water  balance,  more  precise 
thermal  conductivity  variations  must  be  considered  for  accurate 
prediction. 


Fig.  10.  Ratio  of  PCI  molar  flow  rate  to  molar  rate  of  water  generated  in  a  fuel  cell 
operating  at  0.6  V,  1.0  A  cm-2,  353  K  and  101  kPa. 
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6.  Conclusions 

The  relationship  between  thermal  conductivity  and  both 
compression  and  saturation  in  three  representative  fuel  cell  diffu¬ 
sion  media  was  investigated.  Additionally,  the  stress-strain  re¬ 
lationships  of  the  materials  were  measured  using  an  optical 
microscopic  technique.  The  through-plane  thermal  conductivity  of 
diffusion  media  was  found  to  be  highly  dependent  on  the 
compression  and  the  water  content,  and  expressions  relating 
saturation  and  compression  to  thermal  conductivity  were  devel¬ 
oped.  Additionally,  a  theoretical  prediction  of  the  maximum 
through-plane  thermal  conductivity  was  developed  as  a  function  of 
both  compression  and  saturation  to  understand  the  nature  of  liquid 
connectivity  in  saturated  pores,  and  reveals  a  significant  contri¬ 
bution  of  saturation  content  on  conductive  heat  transfer  in  fuel  cell 
DM,  which  should  be  accounted  for  in  multiphase  modeling  effects. 
The  results  of  this  study  should  be  useful  to  improve  predictions  of 
the  internal  temperature  gradients  and  transport  phenomenon 
inside  the  PEFCs. 
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